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Abstract 

We study Klein-Gordon chains with attractive nearest neighbour forces and convex on-site poten- 
tial, and show that there exists a two-parameter family of periodic travelling waves (wave trains) with 
unimodal and even profile functions. Our existence proof is based on a saddle-point problem with 
constraints and exploits the invariance properties of an improvement operator. Finally, we discuss 
the numerical computation of wave trains. 
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1 Introduction 

Chains of coupled particles or oscillators have a broad range of applications in physics, material science, 
and biology, see for instance [DPB93, BK04], and have been studied intensively during the last decades. 
For chains of identical particles coupled by nearest neighbour interactions the law of motion reads 

myj + ivj = ^'{vo+i - Vj) - ^'{yj - Vj-i) - (i) 

Here yj{t) is the displacement of the jth particle at time t, m the particle mass, 7 > a damping 
parameter, and $ and ^' denote the pair and on-site potential, respectively. Examples for such chains 
vifith 7 = arc FPU- like chains with = 0, and Klein-Gordon chains with arbitrary on-site but harmonic 
pair potential. Moreover, in the case ^'(y) — siny one either refers to (1) as the Frcnkcl-Kontorova model 
or the discrete Sine-Gordon equation. 



Major topics in the analysis of atomic chains like (1) are the existence and dynamical properties of 
coherent structures such as travelling waves and breathers, see the review article [IJ05]. A travelling 
wave is a special solution to (1) which satisfies 

y,it)=Yikj~cjt), (2) 
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where k and w denote the wave number and frequency, respectively, and Y is the profile function. 

The existence of traveUing wave solutions to (1) has been investigated by several authors using rather 
different methods. For chains with harmonic <i> we refer to [IKOO] which establishes the existence of small 
amplitude waves by using spatial dynamics and centre manifold reduction, see also [IP06]. Numerical 
simulations are presented in [DEFW93] , and [M JKA02] investigates the existence and stability of standing 
waves by using a continuum approximation for the small amplitude limit. More recently, the existence of 
periodic travelling waves for the Frenkel-Kontorova model was shown in [Kat05] by means of fixed point 
methods, and [BZ06] provides existence results in Sine-Gordon chains with even non-local interactions. 
Close to our approach are [KZ08b, KZ08a], which set the problems also in a variational framework and 
prove the existence of supersonic waves for Sinc-Gordon chains. Finally, a lot of literature addresses the 
existence of travelling waves in FPU-like chains, compare [Pan05, IJ05, Her08] and references therein. 

In this paper we restrict ourselves to Klein-Gordon chains with convex on-site potential ^ and attractive 
nearest neighbour forces, and consider periodic travelling waves, which in turn are called wave trains. 
Due to simple normalisations we can suppose that m = 1 and = ^y'^, so the profile of each wave 
train must solve the nonlinear advance-delay differential equation 

co'-^Y = A,Y-^'{Y), (3) 

where = kj — cot abbreviates the phase, and is a discrete Laplacian with 

iAkY){^) = Yi^ + k) + Y{^ - fc) - 2Y{^). (4) 

Moreover, since (3) is invariant under shifts in (ys-dircction and under the scaling 

Y{ip) -^Y{\lp), k-^X-^k, u-^X-'^LU, 

we can assume that Y is 1-pcriodic with unit cell A = [—1/2, 1/2). 

In order to proof the existence of 1-pcriodic solutions to (3) wc follow a variational approach and charac- 
terise wave trains as solutions to a constraint saddle point problem for the potential energy. In particular, 
the frequency uj turns out to be the square root of the Lagrange multiplier, and cannot be prescribed. 
This is different from the standard variational method which prescribes the frequency and characterises 
travelling waves as stationary points of the action integral, compare for instance [KZOSb, KZ08a]. 

A further key ingredient for our method is an improvement operator, which we introduce below. Due 
to the convexity of 5* this operator possesses nontrivial invariant cones, and this allows to establish the 
existence of wave trains within these cones. More precisely, our main result can be stated as follows. 

Theorem 1. There exists a two-parameter family of solutions to (3) such that the profile function Y is 
1-periodic, unimodal, and even. 

The paper is organised as follows. In §2 we describe our variational setting, and §3 contains the proof 
of the main result. Finally, in §4 we discuss the numerical approximation of wave trains and present some 
simulations. 
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2 Variational setting 



In the remainder of this paper we rely on the foUowing standing assumption. 

Assumption 2. The on-site potential '5 : R ^ R is twice continuously differentiable and uniformly 
convex, in the sense that there exist two constants < m < M < oo such that m < $"(a;) < M for all 
a; G R. Moreover, ^! is always normalised by ^'(0) ~ ^''(0) — 0. 

This assumption in particular implies 

^mx^ < ^'(a;) < ^Mx"^, < M \x\ 

for all x e R, and that 5* has a unique global minimum at a; = 0. We mention that the proofs below 
require m and M to exist only for all bounded sets, so our results remain valid if $" is bounded and 
uniformly positive on each closed interval. 

Spaces and cones of functions We denote by C*^' the space of all functions that are periodic with unit 
cell A and k times continuously differentiable, and equip these spaces with their usual norms. Similarly, 
is the space of all periodic functions which are square integrable on A, and abbreviates the space 
of all functions Y E which have a weak derivative Y' £ L^. Both and are Hilbert spaces with 
scalar products 

{Yi, = f YiMFsMd^, (Fi, y2>Hi = Y^)^, + f ri(^)d^ / Y2{^)d^. 

J A Ja Ja 

Note that (•, is equivalent to the standard scalar product asYt—!- | Y(ip) d(p\ defines a norm for the 
constants. Finally, we denote by Hj the closed subspace of all functions X E with Jj^X{ip) d(p = 0, 
and set 

B^:^{xeK ■■ ill^'ll2<7}- 

To characterise the qualitative properties of wave trains we introduce the cone U of all unimodal and 
even functions on A. that is 

U:= {y e C° : Y{-ip) = Y{ip) and Y{ip) > Y{ip) for aU < < ^ < i} , 

and define a further cone C by 

C:={Y e O ; -Y" eU). 

By construction, we have Y eU \i and only if — + •) eU, and below in Lemma 13 we show that C 
is a subcone of U. 

Examples of wave trains We proceed with some explicit solutions to the wave train equation (3). 
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Remark 3. For k = the wave train equation reduces to the oscillator ODE 

uj^Y" = (5) 



Therefore, there exists a one-parameter family of wave trains with Y € C, which can be parametrised by 



the energy E = WY'"^ + 



Proof. For given E' > wc pick y < such that ^'(t/) = E and define the function Y = Y{t) with 
t G K as solution to Y"(t) = with initial conditions F(0) = y and Y'{0) = 0. According to 

the properties of 5*, compare Assumption 2, this function Y is periodic with period T = T{E), satisfies 
Y{t) = Y{T - t), and is strictly increasing on [0, T/2\. We now define Y &Uhy Y{lp) = Y{Tp + T/2), 
which solves (5) with lo = 1/T. □ 

Remark 4. For the harmonic potential '^{x) = there exists a two-parameter family of wave trains 
with Y G C given by 



N X 9 4sin(fc7r) +c 
y(^) =acos(27r^), u;^ = , (6) 



f 
47r2 

where the wave number k and the amplitude a > are the free parameters. 
Proof. We introduce the Fourier transform {ym),^^^ '^^ that means 



and find (3) to be equivalent to 

^■K^m^u?ym = (2(1 - cos (27rmfc)) + c)y,n, m e Z. 



i(5i 

2 rn 

2sin'^ z we obtain (6). □ 



To solve this we set ym = y-m = \ 5^ni with 5}-,-^ being the Kroncckcr delta, and thanks to 1 — cos2z 



The Lagrangian structure For general and fc 7^ we cannot solve the wave train equation explicitly 
but need more sophisticated arguments to prove the existence of solutions. The starting point for each 
variational approach is the Lagrangian of a wave train 

C{Y) - iv'V{Y) - Vk{Y), r(y) = i / Y\ipf d^, 

with kinetic energy uj'^ViY) and potential energy 

Vk{Y) = \\\VkY\\l+V,,x{Y), Vni{Y) = / ^{Y{ip))dip, 

J A 
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where VkY{f) = Y{f + — y{f — k/2). Notice that Vk is convex due to Assumption 2, and that the 
discrete difference operator Vfc satisfies 



Afe=VfeVfc, VI. = -Vfc (7) 
with Afc as in (4) and * denoting tlic L^-adjoint. 



Our variational method relies on the following main observation. Suppose F e with ij|y||f^i = 7 > 
is a wave train. Then (3) implies that y is a stationary point of Vk under the constraint F < 7, where 
up plays the role of an Lagrange multiplier. To clarify this stationarity condition, we write Y ^ x -\- X 
with a; G M and X G Hg. and restate the wave train equation as 

ui^X" = ^kX -'^\x + X). (8) 



The convexity of Vk now implies that each wave train Y = x + X is a. minimiser for Vk with respect to 
unconstrained variations of x. With respect to variations oi X £ B-y, however, the only minimiser of Vk 
in is the trivial solution X = with multiplier uj^ = 0, and hence we are interested in other types of 
stationary points. Below we show that there exist saddle point solutions Y to (3) which posses a positive 
multiplier cj^ > as they correspond to a maximiser of Vk with respect to variations of X. Moreover, due 
to the properties of the aforementioned improvement operator we can additionally impose the condition 
X e C, and hence we substantiate Theorem 1 as follows. 

Theorem 5. For given 7 > and /c G A there exists a pair [X, a;) G C H S-y x M such that 

Vk(x + X) = max mmVk(x + X). 

The function Y = x + X ^ C is then a wave train, that means it solves (3) for some uj^ > 0. 



3 Proof of the Existence Result 

In this section we always suppose that the parameters 7 > and fc G A are arbitrary but fixed. 

Remark 6. Hq is compactly embedded in and with \\X\\^ < \\X\\^ < \\X'\\2 for all X G Hq. In 
particular, weak convergence in Hj implies strong convergence in both and L^. 

Proof. For given AT G Hq and arbitrary ipo, Lpi the integral representation X{ipi) — X{(po) = J^^^^ X'{lp) 
and Holders inequality imply — ||X'||2 < X{ipi) — X{(po) < ||Ar'j|,, where we used that dip = 1. 
Integrating these estimates with respect to (ySo G A gives < ||Ar'l|2, and \\X\\2 < ||Arj|^ follows 

again from Holders inequality. □ 

Properties of the potential energy We proceed with some elementary properties of Vk- 
Lemma 7. The functional Pfe : ^ M has the following properties: 
1. It is well-defined, nonnegative, continuous, and strictly convex. 



5 



2. It is Gateaux- differentiable with dY'Pk\^] = — A^y + ^'{Y), and its derivative dyVk : ^ is 
a continuous operator. 

3. For all Yi , F2 G we have 

Vk{Y2) - Vk{Y,) >f 11^2 ~ Y4l + {dYV[Y,], Y^ - (9) 

^. y ^ implies Vk{Y) > and dyVklY] ^ 0. 
Proof. The proof of the first two assertions is straight forward. Towards (9) we notice that (7) implies 

WVkYiWl = -{AkYu yi)L2, sIlVfeFaila + ^\\^kYi\\l > -(AfeYi, Y2),_2, 

and hence 

h\^kY2\\l - ^WVkYiWl > -(AfcFi, Y2 - yi)L2. (10) 
Moreover, the convexity inequality for 5* provides 

*(y2)-*(2/i) > f (2/2-yi)' + ^''(yi)(y2-yi), 2/1,^2 eK. 
We set jji = Yi{(p) and integrate this identity with respect to G A to obtain 

VraiY2)~rnl{Yi) > f 11^2 - Yi\\l + {dYVnAYi], Y2 - Y{) . (11) 

The estimate (9) then follows by adding (10) and (11). In particular, exploiting (9) with Y2 ~ Y and 
y = we find VkiY) > for ah Y ^ 0. To complete the proof we suppose that dyVklY] = 0. Then (9) 
with 1^2 = and Yi = y gives 

> -VkiY) = VkiO) - Vk{Y) > f\\Y\\l 

and hence y = 0. □ 

Next we show that for each X we can choose a unique x by minimising the potential energy. 
Lemma 8. There exists a unique and continuous map a; : Hq R such that 

f ^'ixiX)+Xiip))dip = 0, Pk{x{X) + X) = minVkix + X) (12) 

J A 

for allX e Hj. 

Proof. By assumption 2, the function iIj{x) = \E'(a; + X{ip)) dip is well-defined and twice continuously 
differentiable with derivatives ip^^\x) = ^'('^(a; + X{Lpj) d(p for i = 1, 2. In particular, ■(/; is uniformly 
convex with M > il)"{x) > m > 0, and hence there exists a unique minimiser x ~ x{X) with ip'{x) — 0. 
It is straight forward that x{X) depends continuously on X with respect to the strong topology in C°, 
and the compact embedding from Remark 6 implies the continuity with respect to the weak topology in 
HJ. □ 
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In what follows wc consider the reduced potential energy functional 

Pfc : ^ M, Vk{X):^rk{x{X) + X), 

and aim to show that there exists maximisers for Vk in C n ,B-y. To this end we draw the following 
conclusion from Lemma 8. 

Remark 9. The functional Vk is weakly continuous on V\\. 

The improvement operator As a main ingredient for the proof of Theorem 5 we introduce the 
improvement operator "Tfc, 7 as follows. For each X the function X:=7fc_^[X] satisfies 

u^X" = /\kX -^'{x{X)+X), (13) 

where is chosen such that X G dB^. Notice that X is well-defined as (12)i implies that the right hand 
side in (13) vanishes when integrating over A. We further define the integral operator 

{IX)iip):^Jx{ip)d^-J Jx{ip)d^d^, (14) 

AO 

and thanks to (IX)' = X we rewrite (13) as 
TT7 \\T7\\ 

Tk.,[X]:^ ,,2.^NN2^ Z:=dYVk[x{X)+X]=-likX + ^'{x{X)+X). (15) 



Remark 10. The operator T is a well-defined and a compact endomorphism of Hj. In particular, A„ 
Xoo weakly implies TX„ — * lA^oo strongly, andI{X) = implies X = Q. 

Lemma 11. The improvement operator maps B-y \ {0} to dB-y \ {0} and weakly convergent sequences to 
strongly convergent ones. Moreover, it satisfies 

VkiTk,^[X])>rk{X), (16) 

where equality holds if and only if X ~ Tk,^[X]. 

Proof. According to (14) and (15), the function 7fc.^[X] G dB.y is well-defined as long as IZ 7^ 0, which 
holds true if and only if A" 7^ 0, compare Remark 10 and Lemma 7. Moreover, the claimed convergence 
properties are implied by (15) and Remark 10. Now let X £ B^ he fixed, and set X = Tk^j[X] and 
X = x{X), X = x{X). Then (13) reads w^A" = -dyVkix + A), and from (9) we infer that 

Vk{x + X) ~ Vk{x + X) > tj2(-A", + A - A)l2 =a.2^A, a- A)hi. 

This implies (16) due to 27 = ||Aj|^i and (A, A)[^i < ||Aj|^i j|Aj||_|i < 27. Moreover, we have equality in 
(16) if and only if (A, A)^! = ||A||^i = ||A||hi, that means A = A. □ 

The following implication of Lemma 11 is key for our existence proof. 
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Corollary 12. Suppose that the set S C is invariant under the action ofTk^^. Then, each maximiser 
X for Vk in S O By is a wave train with non-vanishing frequency and satisfies ^WXW"^ = 7. 

Proof. Let X be a maximiser for Vk m S O B~^, and recall that X ^ Q according to Lemma 7. By 
assumption, X:=Tk^.y[X] satisfies X G S dB^ as well as Vk{X) > Vk{X), so Lemma ff provides 
:p^^(X)=Pfc(X), and hence X = □ 

Recall that (8) can be viewed as the Euler-Lagrange equation for the optimisation problem Vk{X) 
max with X G B^, where the Lagrange multiplier corresponding to j^X{ip) = vanishes due to the 
choice of x{X). In particular, the fact that each maximiser for Vk in the smaller set 5 n S-y satisfies (8) 
without further multipliers is not clear a priori but provided by the invariance of S under the action of 

Properties of the cones C and U Our next result is rather elementary but provides an important 
building block for our existence result. 

Lemma 13. X G C implies both X eU and — A^X € U. 

Proof. Thanks to Z:= — X" E U the function X is even and there exists i^q £ (0, 1/2) such that 
Z{(p) > for all < (/3 < 930 but Z{lp) < for all (fa < <fi < 1/2. Therefore, X is concave on [—(fa, fo] 
with maximum in 0, and we infer that X'{—(p) > X'{Q) = > X'{lp) for all < <y9 < (po. Similarly, X is 
convex in [ipQ, 1 — ipo] with minimum in 1/2 and thus we have shown X E U. 
Towards the second claim we introduce the averaging operator 

ip+k/2 

and a direct computation shows {AkX {(p)) dip = k Jj^X{ip)dip and that Ak is even provided that X 
is even. Moreover, we have V^X = AkX' = (AkX)' and hence AkX = A^X" for all X. It remains to 
show that AkZ is non-increasing on ip G [0, 1/2] for all Z G U, and to this end we discuss the following 
cases ioT < ip < 1/2: (a) ip - k/2 < < p + k/2 < 1/2, (6) ip - k/2 < < 1/2 < ip + k/2, (c) 
Q<p-k/2<p + k/2< 1/2, [d) < fc/2 < 1/2 < p + k/2. 

In view of Z{p) = Z{—p) we find for case (a) the estimate [AkZ)' {p) = Z(k/2 + (/j) — Z[k/2 — p) < 
0, and Z{p) = Z{1 - p) provides for case (b) that {AkZ)'{p) = Z{1 - k/2 - p) - Z{k/2 - ip) < 0. 
Similarly, for case (c) and (d) we obtain {AkZy{ip) = Z{p + k/2) — Ziip — k/2) < and {AkZ)'{ip) = 
Z{1 — k/2 — p) — Z{p — k/2) < 0, respectively, and the proof is finished. □ 

Corollary 14. X eC implies Tk,j[X] G C. 

Proof. For each X E C Lemma 13 provides —AkX G U and X E U, and from the latter we conclude 
that '^'{x{X) + X) E U as is monotonically increasing. Consequently, the right hand side in (13) is 
contained in — W, and this implies 7fc_-y[X] EC. □ 
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Existence of maximisers The cone C is not closed under weak convergence in Hj, but we can prove 
the foUowing rcsuh. 

Lemma 15. Let X„ d C B-y he any sequence with Xn — > X^o weakly in Hj for some limit X^q 7^ 0. 
Then we have Tk^^[Xn] ^.7 [^00] strongly in Hj and Tk^^lXoc] G C n dB^. 

Proof. For all n G N U {00} let 

Xn:=Tk,^[Xn], Xn:=x{Xn), Zn-^dyV {Xn + Xn) , i^n'-"^' i 

and recall that (15) implies 

LolX'; = -Zn = AkXn - *'(x„ + X„). 

Remark 6 and Lemma 7 provide Z„ Zoo strongly in C'^, and hence also uJn Uoo- Moreover, we have 
'-^00 7^ since X^o ^ implies Zoo 7^ and TZoo ^ 0, compare Remark 10 and Lemma 7. We conclude 
that X!il — > X'^ strongly in C^, thus Xoo G B^, and since U is closed we infer that Z^o G U, and therefore 
XooeC. □ 

Corollary 16. Vk attains its maximum inCHB-y. 

Proof. Suppose that C C H is a maximising sequence for Vk- By weak compactness we 

can extract a (not relabelled) subsequence such that X^ X^o weakly in Hj for some limit X^o and 
Remark 9 implies Vk{Xoo) ~ supT'lcne-,: and hence X^o 7^ 0. Moreover, according to Lemma 15 we have 
^00:=^. 7 [^00] G dB^ n C, and from (16) we infer that X^o is a maximiser. □ 

The combination of Corollary 12 and Corollary 16 gives the proof of Theorem 5. 



4 Approximation of Wave Trains 

By view of the preceding results it seems natural to approximate wave trains by the following abstract 
iteration scheme for fixed points of Tk_^. 

Scheme 17. For given parameters 7 > 0, A; G A and fixed initial value X^ G C H dB^ with ^ Q we 
define sequences 

(x,), CM, K),CR, {X,)^(iCr\dB^ 

by the following recursion: 

1. solve the scalar optimisation problem for Xi = x{Xi), 

2. compute Z^ = -A^Xj + "^'{xi + Xi), 

3. solve Ul' = -Z, forU^eHl, 
4-. compute Lof = \\Ul\\^_2/^/2rf, 
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5. set Xj+i = Ui/ujf. 

From a mathematical point of view the account of this scheme is hmited for the following reasons: (i) 
We have no convergence proof, (n) Due to the lack of uniqueness results it is not clear whether or not 
Corollary 12 covers all fixed points of 7^.^ in C B-y. Nevertheless, suitable discrete variants are easily 
derived and work very well in numerical simulations. 




-0.5 0.0 0.5 -0.5 0.0 0.5 -0.5 0.0 0.5 



Figure 1: Profile X with second derivative and residual for the data from (17). 

In a simple discrete analog to Scheme 17 each profile function Xi is identified with a vector (Xj^)j=i..Ar 
via Xf ~ Xi{—l/2 + j/N), where k is supposed to be a multiple of 1/A^, and the derivative Xl is 
approximated by centred finite differences. Moreover, integrals with respect to ip are replaced by their 
Riemann sums and Xi is computed by a discrete gradient flow for the function x Sj=i jv -^D- 
This numerical approach is illustrated in Figure 1 for the data 

7=10, /s = 0.1, = exp(-a;), (17) 

with TV = 800 and Xo{(p) = ^tt^^ cos {2TTip). Under the iteration the functions Xi converge to a wave 
train, that means a fixed point of 7^, with unimodal and even second derivative. Moreover, numerical 
simulations indicate that the limit profile is independent of the initial profile Xq e C. 



Profile X, k=0.1 Profile X, k=0.3 Profile X, k=0.5 




-0.5 0.0 0.5 -0.5 0.0 0.5 -0.5 0.0 0.5 

Profile V, k=0.1 Profile V, k=0.3 Profile V, k=0.5 




-0.5 0.0 0.5 -0.5 0.0 0.5 -0.5 0.0 0.5 



Figure 2: Profiles X and V for the data from (18). 
Further numerical results are shown in Figure 2 and correspond to 

7 = 50., fc G {0.1, 0.3, 0.5}, ^'"(a;) = l + .T^ (18) 
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where N and Xq are chosen as above. Here we plot additionally the profile function V ~ ~ujX', which 
describes the atomic velocities in a wave train via yj{t) = V{kj — cut), compare (2). 




-0.5 0.0 0.5 -0.5 0.0 0.5 



Figm'c 3: Simulations for the data from (19): Profiles X and V for minimal and maximal value of 7, and 
superposition of six traces of wave trains. 

The simulations in Figure 3 correspond to 

7 e {0.1, 3, 12, 30, 60, 100}, /s = 0.1, ^-"(0;) = exp ( - max{x, of ), (19) 

and illustrate how the wave trains depend on 7. The right picture shows the traces of the wave trains, 
that means the closed curves 

^^{X{^), V{^)), 

whose diameters increase with 7. Surprisingly, we find a nested family of curves: The traces for different 
values of 7 do not intersect, but all traces for 7 < 7 fill out the interior of the trace for 7. This observation 
indicates the existence of a 'hidden structure' for the nonlinear advance-delay differential equation (8). 
In particular, there must be an equivalent planar Hamiltonian system such that the traces coincide with 
the level sets of the Hamiltonian. We cannot prove that the wave trains for fixed k and increasing 7 give 
rise to a nested family of closed curves but mention that a similar phenomenon can be observed for wave 
trains in FPU chains, see [Her08, HR08]. 
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